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Theory of adsorption on a surface with nanolocal defects is proposed. Two efficacy parameters 
of surface modification for nanotechnological purposes are introduced, where the modification is a 
creation of nanolocal artificial defects. The first parameter corresponds to applications where it 
is necessary to increase the concentration of certain particles on the modified surface. And the 
second one corresponds to the pattern transfer with the help of particle self-organization on the 
modified surface. The analytical expressions for both parameters are derived with the help of the 
thermodynamic and the kinetic approaches for two cases: jump diffusion and free motion of adsorbed 
particles over the surface. The possibility of selective adsorption of molecules is shown with the help 
of simulation of the adsorption of acetylene and benzene molecules in the pits on the graphite surface. 
"^h | The process of particle adsorption from the surface into the pit is theoretically studied by molecular 

■ dynamic technique. Some possible nanotechnological applications of adsorption on the surface with 

7-H ' artificial defects are considered: fabrication of sensors for trace molecule detection, separation of 

isomers, and pattern transfer. 
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The progress in the development of nanotechnology gives rise to the possibility to modify the surface at the 
nanometer scale. Various methods of nanometer-scale defects creation using scanning tunneling microscope (STM) 
and atomic force microscope (AFM) were elaborated. In particular, nanometer-scale defects can be produced by 
electric field between the STM needle tip and the surface jij - ||, with the help of nanolocal chemical reactions 
£LT induced by the STM needle tip || - ||, and by direct touching of the STM Rfl or AFM ||] tip. Nanometer-scale 
defects may influence on adsorption of particles (atoms, molecules and clusters) pPf . Surface nanostructures have been 
I ■ fabricated with the help of adsorption on the surface modified by scanning probe anodization |n[]. We believe that 
the further advance in nanotechnology can be achieved when techniques of nanostructures fabrication with the help 
• of adsorption will combine with preliminary nanolocal surface modification. Therefore the theoretical investigations 
t-H ' of the nanolocal processes during adsorption on heterogeneous surface is an actual task at present time. 

Here we present the theory of adsorption on the surface with nanolocal defects, with sizes of the order of the size 
of only one adsorption site (see also JljJl^]). Two efficacy parameters of surface modification by creation of artificial 
defects are introduced. The first efficacy parameter 7 can be useful for such applications where it is important to 
increase the total concentration of particles on the surface. For example, such applications may be fabrication of 
sensors for trace molecule detection and separation of isomers. The second efficacy parameter T can be useful for 
applications where high filling of adsorption sites on defects w 1 is necessary simultaneously with the absence 
of adsorbed particles on the surface free of defects. One of such applications, for example, may be pattern transfer 
through the use of particle self-organization on the modified surface. Two limiting cases of adsorption are considered: 
the case of low temperature when adsorbed particles are localized and the case of high temperature when adsorbed 
O r particles freely move over the surface. Analytical expressions for parameters 7 and T are derived for both cases of 
localized and freely moving particles using equilibrium thermodynamics and the kinetics principles. 
. t-H , The energy of molecular adsorption in pits with various sizes and shapes is investigated. It is shown that this energy 
is determined by the size of the pit. Therefore molecular adsorption on defects may be selective to a type of molecule. 

The process of particle adsorption from the surface into the pit is simulated on using molecular dynamic technique 
on the example of argon atom and graphite surface. The probability of adsorption from the surface into the pit is 
calculated with the help of this simulation. The calculated value of this probability is used to estimate both efficacy 
parameters 7 and T of surface modification for the case of argon atom absorbed on graphite surface with a pit. 

Some possibilities for nanotechnological applications of adsorption on the surface with artificial defects are discussed 
in detail: fabrication of sensors for trace molecule detection, separation of isomers, and pattern transfer as a result of 
particle self-organization on the modified surface. 



*e-mail: lozovik@isan.troitsk.ru, tel: (07095) 334088f , fax: (07095) 3340886 
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II. THEORY OF ADSORPTION ON A SURFACE WITH DEFECTS. 



Let us consider the system containing the particles of one type in three states: in a gas, in adsorbed states on 
defects (adsorption site D) and on areas of the surface which are free of defects (adsorption site S). We restrict our 
consideration to the case of nanolocal defects with only one adsorption site. The process of an adsorption of a particle 
on a surface can be considered as a chemical reaction 

Molecule A + free adsorption site on the surface = 
= adsorption complex A 

The condition for equilibrium of this " chemical reaction" is 

Mg + Ma = Ma (1) 

where /x g , /i a and [i' a are chemical potentials of the particle in the gas, of the adsorbent per one adsorption site and of 
the adsorption complex, respectively. We shall restrict our theoretical consideration to the case of physical adsorption. 
In this case the following approximation can be used pi] 

Ma = Ma + Ms (2) 

where fi s is the chemical potential of the particle adsorbed on site S. We use the analogous approximation for particle 
adsorbed on the defect. Taking into consideration this approximations, the condition for thermodynamic equilibrium 
of the system is 

Ms = Md = Ms (3) 

where fid is the chemical potential of the particle adsorbed on the defect. This case is analogous to adsorption of N± 
immobile particles adsorbed on iV2 adsorption sites. Therefore the chemical potential of the particle adsorbed on the 
defect is (see, for example, Jl4[) 

f M = kT\n Vd (4) 
(1 - v d )qd 

where vs. — Nt/N% is the filling of adsorption sites on defects, qd is the partition function of the particle adsorbed 
on the defect, qd includes internal degrees of freedom, vibrations, and rotations or frustrated rotations relative to the 
surface, and interaction with the surface. The chemical potential of the particle in the ideal gas is 

where qi nt is the partition function of the particle in the gas, qi n t includes internal degrees of freedom and rotations, 
m is the mass of particle, P is the pressure of particles in the gas. In the case of ideal gas the Langmuir isotherm for 
concentration n\ of particles adsorbed on the defects follows from the equation fi g = [id 

q d ex.p(fi /kT)n d P 

m = v d n d = — — (6) 

1 + q d exp(no/kl )P 

where n d is the concentration of defects, and 

, / , f2irmkT\ 3/2 \ 
Mo = -feTln lqr n \ g (_ — J kT J (7) 

For some nanotechnological applications it is necessary to increase the concentration of certain particles adsorbed 
on a surface. Some of these applications, namely, fabrication of sensors for trace molecule detection and separation of 
isomers will be discussed below. We believe that following modification of surface may be useful for such applications: 
On the surface the nanolocal artificial defects are produced so that adsorption of certain particles on defects (sites D) 
causes significant increase in the concentration of these particles in comparison with concentration of these particles on 
a perfect surface (where only sites S are present). Such modification is worthwhile only in the case of low concentration 
of this particles on perfect surface, i.e. when filling vq of sites S is small (vq < 1). By this is meant that the contribution 
of particles adsorbed on sites S into the total concentration of particles adsorbed on a modified surface is negligible. 
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Total concentration nt on modified surface is approximately equal to concentration n\ of particles adsorbed on sites D. 
Let us introduce the surface modification efficacy parameter 7 equal to the ratio between concentrations on modified 
surface n t and perfect surface n 2 

7 = — = (8) 

where n is the concentration of sites S on the perfect surface. It is worthwhile to modify the surface if 7 3> 1. Let 
us obtain the expression for the surface modification efficacy parameter 7 for different cases. The magnitude of this 
parameter is defined by the equation //<j = Ms and therefore it does not explicitly depend on the gas pressure. 

Let us consider the case of low temperature kT <C AE, AE being the height of the potential barrier between two 
adjacent adsorption sites S. In this case particles adsorbed on sites S jump between these sites. The expression for 
the chemical potential of the particle adsorbed on sites S has the same form as for the particle adsorbed on the defect. 
In result the equation \id = /i s leads to the relation 

V A = V A (9) 

(1 - v d )q d (1 - v )qo 

where go is the partition function of the particle adsorbed on site S. Now we substitute in Eq. (@) the equalities 
Vd = ni/nd and vq = n^jriQ and consider the case i/o « 1, when the surface modification is worthwhile for the 
purpose of selective adsorption. In result we get 

— = (10) 

n d - ni n q 

From Eq. @ we obtain 

7 = i Mi(i+^y 1 (ii) 

The partition functions of the particles adsorbed on sites D and S respectively may be represented by the expressions 
@: 

q d = q' d eMEd/kT), (12) 

q a =q' exp(E a /kT), (13) 

where q' d and q' are the partition functions of the adsorbed particles including internal degrees of freedom, vibrations, 
and rotations or frustrated rotations; Ed and Eq are adsorption energies of particles adsorbed on sites D and S 
correspondingly. We get 

- Pnd (14) 



1 + /3n 2 ' 



/?=Aexpfe^) (15) 
n q \ kT J 

Now we consider the case of high temperature kT <C AE when particles adsorbed on sites S freely move over the 
surface. In the case of low concentration of particles adsorbed on perfect surface (i/q <C 1) the interaction between 
adsorbed particles can be disregarded. Therefore we consider the system of particles adsorbed on perfect surface as 
the two-dimensional ideal gas. Then the chemical potential of the particle adsorbed on sites S is given by 

Ms = fcTlnf P^—, ; — ^ (16) 

P \2TrmkTq' s exp{E /kT) J V ; 

where q' s is the partition function of the adsorbed particle, q' s includes internal degrees of freedom, rotations and 
one vibration for constrained motion perpendicular to the surface. From the relation [id = fi s we have the following 
expressions for the parameter 7 
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(17) 



where the parameter /3 for the case considered is 



q' d h 2 ( E d -Ep \ 

Note that in both cases, when adsorbed particles are localized and freely move over the perfect surface, the depen- 
dence of the parameter 7 on the concentrations n 2 and n d has the analogous form (Eqs. (Oh and (J17|)). 



To estimate the value of parameter 7 with the help of Eqs. (14) and ( |17| ) it is necessary to calculate the partition 
functions of the particles. Nevertheless equivalent results can be obtained using kinetics principles. Contrary to our 
derivation of expression for parameter 7 with the help of thermodynamical approach, where we use approximation 
(^) valid only for physical adsorption, the consideration of kinetics below is valid for both physical and chemical 
adsorption. By analogy with the standard premises of the Langmuir and BTE adsorption theories []l7f , we suppose 
that the kinetic coefficients for the adsorption and desorption of the particles are independent of their concentration. 
We suppose also that each site D on a modified surface is surrounded by area with sites S where local concentration 
of particles on this area equal to their concentration on perfect surface n 2 . According to detailed balance for particles 
exchange between D and S sites we have 

k~ d {\-v d )n 2 = k+v d n d (19) 

where kj is a coefficient of the adsorption on site D of the particle which have been just adsorbed on site S, and fcj~ is 
a coefficient of desorption from site D with subsequent adsorption on site S. From Eq. (|l9|) we obtain the expression 
for parameter 7 



n dk d ( 1+ Kt) 
n 2 kj \ kjj 



(20) 



The coefficient is given by the Arrhenius formula fuj| : 

k+ = Q d eM-^E d /kT), (21) 

where Q d is a frequency multiplier equal in order of magnitude to the vibration frequency of the particle adsorbed 
on the defect, and AE d = E d — Eq + AEi is the activation energy for the desorption of the particle from the defect 
onto the surface free of them, E d and Eq being the adsorption energy of the particle adsorbed on sites D and S, 
correspondingly, and AEi is the height of the potential barrier between adsorption sites on the defect and on the 
surface. 

Where kT 3> AE, the particle freely move over the surface, and the coefficient of adsorption on a defect site D 
from an area of surface with sites S is 

k^ = 2k' R <v> s n 2 exp(- AE X / kT) , (22) 

where k! is the probability of the adsorption on the defect site D for the particle moving over the surface just after its 
collision with site D, R is the radius of the defect, and < v > the average velocity component parallel to the surface 
for the particles adsorbed thereon. We define the probability k! as the ratio k' = N a /N c between number N a of events 
of adsorption on site D immediately after the collision and number N c of these collisions. This quantity is analogous 
to the coefficient of attachment of the particles to the surface in the case of an adsorption from a gas on a surface. 
After substitution of Eqs. ( ^l| ) and into Eq. (Effl) one have 

(3n d k'R <v> s { E d -E \ 

-iP= FT ex P ( 23 ) 



' 1 + [3n 2 ' ^ Q d r V kT 

Where kT <C AE, the particles jump between adjacent adsorption sites, and in this case k^ = k'l n u>VQ, where 
k' is the probability of the absorption on the defect immediately after a jump on this site, l n is the number of the 
adsorption sites S adjacent to the defect, uj is the jump frequency defined by the Arrhenius formula 

w = ft exp(--jM , (24) 
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where fio is a frequency multiplier equal in order of magnitude to the vibration frequency of the particles adsorbed 
on the site S. Note that according to simulation fl6f| the majority of jumps are between adjacent adsorption sites. In 
this case 

Pn d p_ k'l n n ^ / E d -E \ 
1 + f3n 2 ' n a tt d \ kT ) 

Expressions Jl4|), (|T^), (pg]), ( p^ ) can be extended to the case of some different kinds of defects. Disregarding the 
contribution of particles adsorbed on sites S to the total concentration n t on of particles adsorbed modified surface 
we get 



n, 



nu (26) 



nii is the concentration of particles adsorbed on the i-th kind of defects. In this case the surface modification efficacy 
parameter 7 for the total surface is 



7 



n t /n 2 = ^n li /n 2 = 7i (27) 



where 7, is the surface modification efficacy parameter for the i-th kind of defects. 

For some nanotechnological applications a high filling v d ~ 1 of sites D and absence of adsorbed particles on sites 
S are necessary. For example, let us consider a pattern drown by the particles adsorbed on sites D. To distinguish 
such pattern two conditions should be fulfilled. Firstly, a high filling v d w 1 of sites D corresponds to a large value of 
parameter ot\ — 1/1 — v d . Secondly, for better pattern contrast the filling v d of sites D should be significantly greater 
than the filling vq of sites S, that corresponds to a large value of parameter a 2 — v d /vr>- Let us introduce the second 
surface modification efficacy parameter T which is equal to product of parameters a\ and a 2 

T = ai a 2 = —pt (28) 

fo(l - v d ) 

The analytical expressions for parameter T are derived with the help of both thermodynamical and kinetic ap- 
proaches. Under the thermodynamical consideration we use the equation /i d = fj, s analogous to the derivation of 
expressions for parameter 7. In the case of jump diffusion of particles adsorbed on sites S we get 

r q' d E d - E 

T = -fexp — — — (29) 
q' kT 

When the particles adsorbed on sites S freely move over the surface we get the following expression for parameter 

r 

n h 2 q' d E d - E , . 

r = o exp — ur — ( 30 ) 

Zirmkl q kl 

The Eqs. (|9|) and (|o|) are not convenient for practical use (analogous to Eqs. ( |l~4|) and (|l7|)) because it is necessary 
to calculate the partition functions. Therefore we have derived the equivalent expressions with the help of kinetic 
consideration. In the case of jump diffusion of particles adsorbed on sites S we get 

k'l n n E d -E Q 
T = — e * P ^f- (31) 

In the case when particles adsorbed on sites S freely move over the surface we obtain the following expression 

k'R < v > s n E d - E Q 

r = n d e ^^f- (32) 

Note, that nanotechnological procedures described by parameters 7 and T are applied only to the case of low 
concentration of adsorbed particles on sites S (v <C 1). Therefore the neglect of interparticle interaction is adequate 
for our consideration. However, Eqs. (p3|), (p5|), (^Tj) and ( |32| ) derived in the framework of the kinetic approach may 
be easily extended to the case where this interaction is important (analogous to homogeneous adsorption, see, for 
example, ]T^|). To take into account an interparticle interaction it is sufficient to replace E d by E d + w d N d and E by 
Eq + woNq, where w d and wq are interaction energies between neighbor particles, and N d and Nq are average numbers 
of nearest neighbors for particles adsorbed on sites D and S, respectively. 



5 



III. ADSORPTION ENERGY OF MOLECULES IN THE PITS. 



The adsorption of particles on an unperfect surface for the simplest cases had been studied previously: in the 
cylindrical pore Jl9| , in the split-shaped pore [ pp[ , at the intersection of two perpendicular steps forming a reentrant 
corner pp| ]. Nevertheless these studies dealt only with rather large surface defects. The case of nanolocal defects with 
sizes of order of the size of only one adsorption site was not considered. 

Theory presented above does not depend on the nature of defects. They may be structural defects or chemical one 
and so on. We consider molecule adsorption on a single type of structural defects, namely, pits made by removing 
tens of atoms from the surface. Our interest in this type of defects is attracted by following reasons. Firstly, according 
to expressions (fL4|), (|l7|), (|23|), ( p5| ) surface modification efficacy parameter 7 ~ exp((i? p — E$)/kT). Therefore the 
modification of surface by presence of the pits with certain size can cause the increase of concentration of certain kind 
of molecules. Some nanotechnological applications based on this effect are discussed in Sec. 4 (fabrication of sensors 
for trace molecule detection, separation of isomers etc.). Secondly, because of the recent nanotechnological advances 
the various techniques to produce pits on the surface had been elaborated. For example, it is possible to fabricate 
the pits only several nanometers in diameter and even to remove a single atom from the surface structure || 
with the help of field desorption by STM tip. Therefore a theoretical study of particle adsorption in a pit is an actual 
problem now. 

Here we have calculated the adsorption energy E p for the argon atom and benzene and acetylene molecules adsorbed 
on graphite surface in pits that are different in size but can accommodate only one particle. We have choose a graphite 
for our study because the adsorption of argon (see, for example, ^J) and various molecules (see, for example, f2^] 
p4|) on perfect graphite surface have been carefully theoretically studied in a set of papers. Besides, the pits on the 
graphite surface were produced by nanolocal chemical reactions near STM tip ^ and by bombardment of cations 
with subsequent etching p5{ ]. 

This is a nonspecific adsorption case where the adsorbate-adsorbent interaction is of the Van der Waals type jl7) . 
We have described the interaction between the argon atom and the carbon atoms in terms of the Lennard- Jones 
potential 

U = Ae[{a/rr-{a/rf] (33) 

with a cutoff radius of r c — 3.2, where the potential parameters have been taken to be a = 3.12 A and e = 54.4 K 
[2l|, a = 3.82A and e = 31.6 K for the C-C interaction |||, and a = 3.37 A and e = 21.7 K for the H-C interaction 
[ 26 1 . The atoms of graphite have been assumed to be fixed at the lattice sites, three graphite layers with 288 atoms 
in each layer being taken into account. The interatomic bonds in the molecules have been taken to be rigid, and the 
interatomic bond angles fixed. The bond length values used (1.4 A for the C-C bond and 1.08 A for the H-C one) 
have been borrowed from |2q ] . 

We have considered pits with near circular shape that are created by the removal of atoms located inside circles of 
radii R\ and R2 for first and second graphite layers correspondingly. The centers or these circles locate on one vertical 
line I. All possible pits with < R2 < Ri < 5. 06 A for argon atom and with < R2 < Ri < 6. 80 A for acetylene 
and benzene molecules have been investigated. Several different positions of line / relative atoms of graphite were 
considered. Namely, the line I passes through: 1) an atom of upper graphite layer under which an atom of the second 
layer locates (type A); 2) an atom of upper graphite layer under which an atom of the second layer is absent (type B); 
3) the center of a bond between two atoms of upper graphite layer (type C); 4) the center of hexagon formed by atoms 
of upper graphite layer (type D). The optimal pits with one and two graphite layer deep, where the adsorption energy 
is a maximum E™ ax are found for all three adsorbates. For all these cases the energies of adsorption in optimal 
pits are maximum when the center of a particle coincide with the center of the pit. Several pits of oblong shape 
with the size close to the size of the optimal pit with near circular shape are also examined for acetylene molecule. 
Nevertheless the adsorption energies in these oblong pits are less than adsorption energy in the optimal pit with near 
circular shape. These adsorption energies E™ ax and characteristics of optimal pits are listed in Table 1 for all three 
adsorbates. 

The calculations show that the adsorption energy in an optimal pit for small particles (argon atom and acetylene 
molecule) is nearly twice as large as the adsorption energy on perfect surface. The adsorption energy in the optimal 
pit with two graphite layer deep is only 7 % greater than the adsorption energy in the optimal pit with two graphite 
layer deep. The further increase of pit depth does not cause the adsorption energy to increase. Therefore a pit a mere 
1-2 atomic layers deep may be sufficient to increase the absorption energy to the utmost, which is fairly convenient 
for nanotechnological applications. 

A particle adsorbed in a pit interacts with walls and bottom of the pit. And a particle adsorbed in a capillary 
interacts only with walls of capillary. Therefore the adsorption energy of a particle in a pit is greater than in a 
capillary with the same radius. Previously the adsorption energy of a particle in a cylindrical capillaries with various 
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radii has been calculated analytically [|19|. In |19| the interaction between the particle and the walls of the capillary 
has also been described in terms of the Lennard- Jones potential. However, the capillary walls in this simplified model 
have been considered continuous. According to calculations in the framework of this model the ratio eo between the 
maximal adsorption energy E™ ax of a particle in the capillary and the adsorption energy on the perfect surface is 
6q = E™ ax /E s ~ 3. This value of eo is 1.5 times as great as value calculated by us for adsorption in a pit eo = 2.06. 
Therefore the consideration of the location of all the atoms of surface is necessary for exact calculation of adsorption 
energy in a pit or a pore. For benzene molecule adsorbed in an optimal pit, the simplified model [ fl9| overestimates the 
maximal increase in the adsorption energy even more because atoms of adsorbed molecule interact weaker with pit 
wall sections distant from them than with the pit bottom. The value of the energy ratio eo calculated here for argon 
atom adsorption in a pit on graphite surface is close to the one calculated for argon atom adsorption on microporous 
magnesium oxide with taking into account location of all the atoms |2(|. According to the calculation |2(J E™ ax /E s 
~ 1.8 for the split-shaped pore and at the intersection of two perpendicular steps forming a reentrant corner. 

Fig. 1 demonstrates that the energy of adsorption in a pit is determined by the size of the pit. This figure also 
shows that a pit too small to accommodate a molecule even reduces its adsorption energy. This two facts open up 
new interesting possibilities for nanotcchnological applications which have been discussed below. An additional point 
to emphasize is that the "flat" benzene molecule placed in a pit of optimal size but with a "unflat" bottom has a 
lower adsorption energy than in the case of unmodified surface. 

IV. SIMULATION OF ADSORPTION OF PARTICLE IN PITS. 

The parameters 7 and T have been estimated for the case of argon atom adsorption on graphite surface with 
pits. For this estimation it is necessary to calculate the oscillation frequency fl p of atom adsorbed in a pit and the 
probability k' of the absorption in the pit for the particle moving over the surface just after its collision with the pit. 



These quantities were obtained by molecular dynamics simulations. According to Eqs. (23) and ( p2| ) parameters 7 
and r increase with the adsorption energies E p in a pit. Therefore two pits were used for simulation: the optimal 
pit with maximum adsorption energy in it and the pit with maximum adsorption energy among the investigated pits 
with only one graphite layer depth (see Table 1). The adsorption energy in this pit is only 8% less than adsorption 
energy in optimal pit. 

A system consisting of three graphite layers with 288 atoms in each layer was used in simulation. Along X and Y 
axes parallel to the graphite surface we impose periodic boundary conditions on the system. The size of the simulation 
cell in these directions was 27.08 x 23.43A The interaction between the argon atom and carbon atoms was represented 
by Lennard- Jones potential (|33|). The atoms of second and third graphite layers were fixed at equilibrium positions. 
The interaction between carbon atoms of first graphite layer was described by modified Born potential 

a N 1 \ a N 

u=^y: r j ~ ro " )r " ) 2 +~e ('« - 04) 

where are distances between carbon atoms, roy are distances between carbon atoms at equilibrium positions, a 
and (3 are force constants. We take a = 505.1 N/m and (3 — 84.4 N/m p7fl . The Born potential represents the 
expansion of the interaction energy between carbon atoms in terms of Ar^ = — . Therefore this potential is 
adequate only for small values of Ar^. We simulated the system at temperatures that approximately 20 times less 
than the temperature of graphite melting. At such temperatures the displacements of graphite atoms from equilibrium 
positions and, consequently, values of Ar^ are small. Therefore we consider that Born potential is acceptable for our 
simulation. 

The equations of motion were integrated using the leap frog algorithm. The integration step used was r = 2. • 10~ 15 
s. Initially the system of graphite atoms has been brought to the equilibrium during 5 TO 3 steps (about 300 oscillations 
of atoms in graphite) in canonical ensemble and 5 • 10 3 steps in microcanonical ensemble. Further simulations were 
performed in microcanonical ensemble. The total energy of the system with the temperature 185 K was conserved to 
within 1 % and average fluctuations of temperature were within 3 %. 

The estimation of oscillation frequencies fid of atom adsorbed in the pits has been performed for the system with 
the temperature 165 K. Initially we took the system without a pit with equilibrium coordinates and velocities of 
carbon atoms. Then the pit arises near the center of simulation cell and an atom was placed in the pit at the position 
that corresponds to the maximum adsorption energy. The new system with the pit and the atom adsorbed in this pit 
has been brought to the equilibrium during 2 • 10 3 steps. Then frequencies fid were calculated during 6 • 10 3 steps. The 
obtained values averaged over 40 modeling experiments are fldi = 5.3 ■ 10 11 s _1 and fid — 4.7 • 10 11 s _1 for optimal 
pits with one layer and two layers deep. 
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The size of simulation cell used is not sufficient for argon atom appeared in any place of the cell to come to 
equilibrium with the surface before a collision with the pit. To overcome this difficulty the simulation of the process 
of the atom adsorption from the surface into the pit was carried out using the following procedure. An argon atom is 
placed in the center of simulation cell without a pit at the position that corresponds to the minimum of interaction 
energy with the surface. Then this system comes to equilibrium until the argon atom reaches the boundary of the 
cell. At the instant the argon atom crosses the boundary the pit arises in the center of the cell. The appearance of 
the pit was investigated for the system at the temperature 185 K. The change in the energy of interaction between 
argon atom and graphite atoms caused by the appearance of the pit is within 0.03%. Therefore we consider that the 
appearance of the pit does not essentially perturb the motion of argon atom and equilibrium between argon atom and 
surface. During the time between the appearance of the pit and the argon atom collision with the pit graphite atoms 
make about 100 oscillations that is sufficient for system to come to equilibrium p8| ]. 

On the average an argon atom makes 12 jumps before collision with the pit. Average change in the kinetic energy 
of argon atom in result of collision with the surface is < AE >= 17 ± 1.2%, AE, = 2\E t - • 100%/ (Ei + £?i-i), 

where £7j_i and Ei are average through the time of jump kinetic energies of argon atom between i-l-th and i-th, 
and i-th and z+l-th between collisions with graphite surface respectively. The value < AE > has been obtained by 
averaging over 20 argon atoms, each making 10 jumps along the surface. We believe that observed energy exchange 
between argon atom and the surface is sufficient to argon atom to come to equilibrium with the surface. 

An argon atom was considered to collide with the pit when the distance between it and the center of the pit was 
less than the distance between the center of the pit and a adjacent adsorption site S. An argon atom was considered 
to be adsorbed in the pit when it makes two oscillations in the pit. The simulation with one atom was performed 
until the atom reflects from the pit, adsorbs into the pit or evaporates from the surface. To exclude atoms with low 
velocities the time of experiment was limited to 2 • 10 4 integration steps. The results of simulations are averaged over 
all modelling experiments for given temperature of surface. These results and calculation of the probability k! of the 
adsorption into the pit and the efficacy parameters are presented in Table 2. The probability k' of the adsorption 
into the pit slightly increases with temperature. We offer the following explanation for this increase. When a particle 
occurs on the surface near a pit edge it interacts with a less number of atoms of the surface than a particle located 
on perfect surface. Thus the pit is surrounded by an energetical barrier. The magnitude of this barrier is found to be 
about 300 K in simulated system. Therefore the fraction of particles capable to overcome the barrier increases with 
the increase of temperature for the considered temperature range 120-280 K. 

According to our estimations for considered case of argon atom and optimal pit in a graphite surface at temperature 



200 K we get pn < 1 for v < 1. Then Eq. (|23|) is simplified 

7 = find (35) 

so that the parameter 7 does not depend on the concentration of atoms in the gas. The Eq. ( |35| ) have been used 
to estimate the parameter 7 for considered case. Let us express the concentration of pits in terms of fraction of the 
surface occupied by them, i.e., = a/irR 2 . The parameter 7 is calculated for a = 0.1. The calculated parameters 7 
and r are presented in Table 2. The magnitudes of these parameters are sufficient for nanotechnological applications 
discussed below. We believe that this conclusion may be also obtained for various other pairs adsorbent-adsorbate. 

Graphite is a very homogeneous adsorbent, i.e., the potential barriers AE between the adjacent adsorption sites 
on it arc low. For the argon atom, AE/k « 32 K ]29|] . At low temperatures (kT < AE) the submonolayer [vq <C 1) 
adsorption of argon on a graphite surface is impossible for any experimentally obtainable concentrations of atoms in 
the gas. Therefore the case of jump diffusion of adsorbed particles is not considered for this pair adsorbent-adsorbate. 

V. NANOTECHNOLOGY APPLICATIONS. 

Here we describe some possible nanotechnology applications for the adsorption on a surface with nanometer-scale 



artificial defects (see also 1 30 ) : pattern transfer as a result of self-organization of particles deposited on such surface, 
and using of selective adsorption for sensors fabrication and isomer separation. 

Pattern transfer, he elaboration of methods of pattern transfer in nanometer scale, i.e. fabrication of nanos- 
tructures on the surface according to a given scheme is very important for some nanotechnology applications (for 
example, superdense recording of information, fabrication of periodic arrays of quantum dots and quantum wires, 
one-electron devices etc.). In principle, surface nanostructures had been produced with the help of STM or AFM tip 
by the transfer of single atoms | ]3l| , p2[ , clusters |33||34||, and nanoparticles |Q. Nevertheless, until now these methods 
have the productivity insufficient for nanotechnology purposes. 

The other set of methods for surface nanostructures fabrication is based on self-organization of deposited particles: 
a) the formation of islands in the result of nucleation in two-dimensional film (see, for example, Jl(|); b) self-organizing 
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ordering in epitaxial layer with facets formation (see, for example, |3q]); c) the formation of periodical surface structure 
by laser radiation [ p7|j3q| . The method a) allows to control only average distance between islands. The methods b) 
and c) are suitable only for periodical nanostructures fabrication with periods by chemical composition of nanoobjects 
and wave length of laser radiation, correspondingly. The productivity of these methods is considerably greater than 
transfer of single particles by STM or AFM tip. Nevertheless they are not suitable for pattern transfer. 

A particle adsorbed on defect may play the role of nuclcation center for island formation from the particles deposited 
on the surface. For example, the gold islands on the graphite surface irradiated with ions form only around defects 
fiof . We propose to deposit particles on the surface with defects located with the help of STM or AFM according to 
given scheme. In this case the formed islands would merge into planned surface nanostructure. 

Recently the diode on one molecule was fabricated |59) . We proposed to place such molecules in electronic schemes 
of nanometer scale with the help of adsorption on specially created defects. We believe that by these means it is 
possible to control not only the position but also the orientation of the molecule. 

Sensors for molecular detection. Methods for detecting traces of certain molecules in a gas where their 
concentration is low are of great importance in modern technology and find application, for example, to ensure safety 
in chemical industries and monitor environmental pollution. The further development of such methods, particularly 
is therefore a high-priority task. A weak spot in methods used is the collection of molecules from air. The known 
method of improving sensitivity [fjfj|| involves the nonselective accumulation of the trace molecules of interest on a 
cooled substrate and their subsequent pulsed laser desorption. It therefore seems extremely tempting to develop a 
method for selective adsorption of molecules on a cold surface. 

Calculations performed show that the adsorption energy in a pit is determined by the size of the pit. We propose 
to modify the surface of sensor by presence of the pits with maximum adsorption energy for certain kind of molecules. 
According to the theory presented the concentrations of these molecules on the surface of sensor will increase. Therefore 
it would enable one to add one more selective step to the existing two detection steps — the optical and the mass- 
spectrometric ones. 

Isomer separation. Pits can be arranged on a surface so closely that the distances between them have the same 
order as their size. The molecular adsorption energy on such a surface will be higher for molecules that fit in the 
pits, compared to that on the unmodified surface, and lower for those which fail to fit in the pits. A surface modified 
in this fashion could be used to effect the selective adsorption of a particular molecular species from a mixture of 
different molecules. What is more, a surface can be modified to have pits capable of accommodating only one of 
several isomeric molecular species. We suggest using such a modified surface to detect or separate molecules differing 
in shape only (e.g., to isolate linear or cyclic isomers from their mixture and separate fullerene isomers and left- and 
right-handed molecules) . 

In summary, the theory of adsorption on the surface with nanolocal defects is developed. Two efficacy parameters 
of surface modification by creation of artificial defects are introduced for different nanotechnological applications. 
The estimations with the help of molecular dynamics simulations on the example of argon atom adsorption in the 
pit on the graphite surface show that magnitudes of these parameters are sufficient for possible nanotechnological 
applications. 
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Figure captions. 

Fig.l. The dependence of the energy E p of adsorption in a pit with one graphite layer deep on the number Ni of 
atoms removed when making this pit; a) acetylene molecule, b) benzene molecule. Solid squares corresponds to pits 
of type A and open circles correspond to pits of type C. Points corresponds to pits of types B and D are not shown 
because they couincide with points corresponds to pits of types A and C. 
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Table 1. Characteristics of optimal pits: the type of the pit I (see text); the number of atoms removed when 
making this pit from the first and second graphite layers, N\ and A/2, respectively. Adsorption energies in the optimal 
pit E™ ax and in the optimal pit relative to that on the unmodified surface AE = E™ ax — E$ (in Kelvin degrees); the 
energy ratios e = E™ ax / E and ei = E™ ax /E™[ ax , where E™{ ax is the maximum adsorption energy among the pits 
with one graphite layer deep. 
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Table 2. Results of simulation of adsorption of atom from the surface into the pit and calculation of surface 
modification parameters, h is the depth of pit in graphite layers; T g is the temperature of graphite surface; AT C is 
the number of atoms collided with the pit; A a is the number of atoms adsorbed into the pit; k' is the probability 
of the absorption into the pit, 7 and T are surface modification parameters. The accuracy of calculation of k' and 
parameters 7 and T is connected with statistical error of detected events of atom adsorption into the pit. 
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